SYSTEM AND METHOD OF DATA REDUCTION FOR IMPROVED EXPONENTIAL 
DECAY MEASUREMENTS 

RELATED APPLICATIONS 

5 This application incorporates by reference U.S. patent 

number 5,52 8,040 to Lehmann in its entirety to the extent it 
is not inconsistent with the application presented herein. 

FIELD OF THE INVENTION 

10 The present invention is directed toward a data 

processing system and method where a signal is contaminated 
by noise of a frequency much less than those that dominate 
the expected signal itself. More specifically, the present 
invention is directed toward a system and method for 

15 optimized data processing in ring-down spectroscopy. 

BACKGROUND OF THE INVENTION 

Ring- down spectroscopy determines the concentration of 
an absorbing species within an optical resonator by the 

20 measurement of a single parameter, the exponential decay 
rate of radiation in that resonator. Invariably, regardless 
of the experimental details of the particular scheme that is 
used, the measurement consists of detecting the radiation as 
it decays with an optical square law detector, amplifying 

25 the detected signal, and then processing the resulting 
record to determine the decay rate. 

Because of the presence of fundamental broadband 
(white) noise in measurements, whether dominantly from the 
30 detector, the amplifier or the light itself, fluctuation and 
error occur in the measurement and the decay rate must be 
determined by one of a number of possible statistical 
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regression (curve fitting) techniques. Often, a single 
measurement is insufficient for the desired precision so 
instead the measurement is repeated many times to produce an 
ensemble of records from which a more precise averaged value 
of the desired parameter can be obtained, provided there is 
no drift in the measured quantity. In the averaging process 
the improvement in precision is proportional to the square 
root of the number of records. In a practical instrument, 
the speed of measurement is important and the processing 
that achieves the desired precision with the fewest records 
(i.e. the shortest time) represents the preferred process. 

For instruments operating in various real -world 
industrial environments there may often be sources of noise 
beyond those noted. For example, the ubiquitous sixty-Hertz 
line frequency from transformers, imbalanced grounds, and 
close proximity heavy duty machinery, can be picked-up and 
produce noise spectral density well in excess of the level 
of fundamental white noise. Careful design of the 

electronics package has been shown to reduce, but rarely 
eliminate, the effects of external noise sources. Numerical 
simulations have established that sixty-Hertz noise signals 
with a root mean square ("RMS") amplitude comparable to that 
of the broadband noise, but unsynchronized with respect to 
the initiation of the ring-down signal, can more than 
quintuple the uncertainty in the measurement, resulting in a 
nearly twenty- five fold increase in the number of records 
needed to achieve the desired precision. 

Therefore, it is an object of the present invention to 
provide a system and method that greatly reduce the effect 
of such noise, reducing the effect for the case sited to a 
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mere one and a half fold increase in the number of records 
needed to achieve the desired precision. 

Additionally, it is another object of the present 
5 invention to achieve a significant improvement in the 
precision of filtering any excess noise components that have 
a period that is greater than four times the length of a 
record. 

10 SUMMARY OF THE INVENTION 

The present invention overcomes the limitations of the 
above prior art techniques for measuring exponential decays 
in the presence of excess low frequency noise. 

15 The present invention is an efficient addition to 

existing systems and methods for data processing of ring- 
down signals that correct for errors introduced by low 
frequency noise that may be present, e.g., due to line 
frequency pick-up. 

20 

In the present invention, the required length of a 
record is determined by a number of factors, including the 
empty resonator decay time and the ratio of the signal 
amplitude to the RMS noise background. Typically a record 
25 length is several times longer than the time for the 
exponential signal to decay to within the RMS level of the 
background white noise. 

The best fit parameters of a function that is nonlinear 
30 cannot be directly calculated in general, forcing the use of 
iteration techniques, that require initial guesses for the 
parameters, that can consume large amounts of computer time, 
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and ultimately can fail to find a result. Such a limitation 
of the prior art is described in "Numerical Recipes," 
William H. Press, et al, page 521, Cambridge University 
Press, 1986, New York, NY. Fortunately, in the case of an 
5 exponential decay to a known background level, it is 
possible to transform the function by a logarithmic 
transformation into a linear one that easily can be fit. 
This technique is described in "Data Reduction and Error 
Analysis For The Physical Sciences," Philip R. Bevington and 

10 D. Keith Robinson, page 134, Second Edition, McGraw-Hill, 
Inc. New York, NY. This method can be applied to the 
digitized data after the subtraction of any DC level that is 
present from the electronic processing. That level may be 
obtained at the beginning of the measurement process by 

15 taking a measurement with no ring-down signal or preferably 
on each record to eliminate possible drift by measuring the 
DC offset far into the exponential decay where the signal 
has fallen significantly below the background noise 
amplitude. The variance in that level is used as an estimate 

20 for the variance of the data points and in the subsequent 
weighting of the linearized fit, as shown in "Pulsed, 
single-mode cavity ring-down spectroscopy," Roger D. van 
Zee, et al, Applied Optics, 38, 3951(1999). 

25 Once the data has been conditioned by removal of the DC 

level, an additional conditioning is necessary before the 
logarithmic transformation can be made. Data points that 
are of zero value or negative have no real logarithm and are 
thus, removed. These data points first appear once the 

30 exponential decaying signal is comparable with the noise and 
define the data cutoff point for fitting. In addition, a 
few of the initial data points are also removed because of 
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transients in the radiation switching. The estimates for 
the parameters are then obtained from the weighted linear 
least squares fit to the logarithm of the conditioned data. 
Here, where y = A exp{-rt}, the amplitude, A, and the decay 
5 rate, r, are the parameters. The variables are y, the 
voltage above background and t, the time after initiation of 
the ring-down. 

When only white noise is present, the resulting 
10 estimates for the parameters and their variances follow the 
theoretical statistical predictions. However when a 

relatively low level of excess noise is present at low 
frequencies, there is a significant increase in the variance 
of the parameters. Further, the standard statistical 
15 prediction for the uncertainty of the fitted decay rate 
substantially underestimates the true uncertainty, as found, 
for example, by looking at the variance of the decay rates 
extracted from an ensemble of identical (except for noise) 
decay events. 

20 

In general, one possibility to compensate for excess 60 
Hz noise alone is to include in the fitting routine an 
additional sinusoidal function that is a linear composite of 
a sine and a cosine function. Unfortunately, this addition 
25 to the overall equations would make them fully nonlinear and 
thus drastically increase the computational resources 
required to obtain the desired solution. 

An efficient simple solution of the present invention 
30 is available by recognizing that the effect of a low 
frequency sinusoidal signal to the first order is to produce 
an overall slope to the record. The present invention fits 
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the tail of the record to a straight line. The line is 
extrapolated to the entire record and then subtracted from 
the initial data points before the logarithmic 
transformation is performed. As long as the record is 
5 smaller than a quarter period of the sinusoid, a significant 
improvement is obtained. This makes the system and method 
more robust since it will correct for low frequency noise 
regardless of the exact frequency of that noise, or even if 
it is made up of several frequency components, provided the 
10 period of the highest frequency is greater than four times 
the record length. 

BRIEF DESCRIPTION OF THE DRAWINGS 

15 Exemplary embodiments of the invention are described 

below with reference to the following accompanying drawings. 

Figure 1 is a schematic diagram of a sample apparatus 
suitable for implementing the methodologies and systems of 
20 the present invention. 

Figure 2 is a flow diagram illustrating certain 
methodical aspects of the present invention. 

25 Figure 3 is a graph depicting a data model of exemplary 

ring-down signal sequence and noise. 

Figure 4 is a graph depicting a data model of exemplary 
ring-down signal sequence and noise in component parts. 

30 

Figure 5 is an expanded view of the first ring-down 
signal in Figure 3 . 
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Figure 6 is an expanded view of the second ring-down 
signal in Figure 3 . 



5 Figure 7 is an expanded view of the third ring-down 

signal in Figure 3 . 

Figure 8 is an expanded view of the first ring-down 
signal in Figure 3 . 

10 

Figure 9 is an expanded view of the first ring-down 
signal in Figure 5 after the DC baseline has been 
subtracted . 

15 Figure 10 is an expanded view of the first ring-down 

signal in Figure 5. 

Figure 11 is an expanded view of the first ring-down 
signal in Figure 5 showing a fit line that will be 
20 subtracted from the original data. 

Figure 12 is a graph of a ring-down signal of Figure 11 
in an expanded view. 

25 Figure 13 is a graph on a logarithmic scale of the data 

in Figure 12 . 

DETAILED DESCRIPTION OF THE PRESENT INVENTION 

The present system and method has been tested with 
30 experimental data and with data generated from a computer 
model. In each case, a significant improvement was 
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measured; with the computer model the improvement was 
quantified . 

Appendix A is provided as the source code used in 
5 generating the computer data. A model of a ring-down signal 
with noise is created to mimic experimental data obtained 
from a photodector and preamplifier applied to an analog to 
digital converter. It is then fit with the new and improved 
algorithm that subtracts a line that is fit to the 
10 background before applying the standard algorithm and 
compares that fit with one that simply uses a constant for 
baseline subtraction. 

For the computer program, exponential decays of tau = 
15 8 0 microseconds were used. A peak signal was chosen. 
Uncorrelated Gaussian noise was added. And a dc offset and 
sixty Hz noise were added. As in real experiments the data 
is sampled every microsecond. 

20 The following section illustrates results from the 

source code written for the algorithm of the present 
invention. This source code, given in Appendix A, is written 
in the Fortran 77 Programming Language. The specific type 
of high level language used to implement the algorithm is 

25 not particularly relevant to the present invention, since it 
will be apparent that almost any type of high level language 
such as Microsoft's VISUAL BASIC™, VISUAL C or C++ could be 
used instead. The following brief example of the application 
of the algorithm to synthetically generated data 

30 demonstrates the improvement in the decay time variance with 
the present invention. 
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Tests were conducted using the system and method of the 
present invention. The test results were computed using 500 
independent ring -downs with the RMS of the 60 Hz equal to 
the RMS of the Gaussian white noise. These results are 
5 provided as Appendix B. As indicated on the bottom of the 
table in Appendix B, the squared ratio of the ensemble sigma 
tau is 19.77. That ratio shows that to achieve a given 
precision in tau, in the prior art, one needs 19.77 times as 
much averaging using the standard algorithm compared with 

10 the new system and method of the present invention. The run 
of 500 independent ring-downs was repeated ten (10) times to 
develop a statistic on the improvement for this set of 
parameters and the result was that the average improvement 
of the ensemble sigma tau was by a factor of 18.98 with a 

15 standard deviation of 1.42. As expected, if the 60 Hz noise 
signal gets smaller, the improvement also gets smaller. 
However, the significance of the result is that the time 
taken to make a measurement of a given precision is reduced 
by that ratio. To be most effective, the low frequency 

20 noise that the algorithm corrects for should be less than or 
equal to 0.42/Tr, where Tr is the length of the record in 
seconds . 

Chosen values for the sample are as follows: amp = 
25 800.00; 60 Hz-RMS = 3.50 offset = 10.00; tau = 80.00; and 
gauss RMS = 3.5. These values were typical of those measured 
in the laboratory, where the amplitudes and voltage levels 
are measured in millivolts, and time in microseconds. 

30 The output of the source code produces the table in 

Appendix B. The resulting values follow: average tau = 
80.13, sigma tau = 0.8777, average reduced (chi) 2 = 1.099, 
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average amp = 799.27, sigma amp = 1.0028, and the squared 
ratio of sigma tau from the traditional and new algorithm = 
19.78 . 

5 An example of a cavity ring-down spectroscopy (CRDS) 

system is shown in Figure 1, as taught in U.S. patent number 
5,528,040 by Lehmann, to provide a better understanding of 
the present invention. The apparatus of Figure 1 can be 
used in conjunction with the present invention, although 

10 other variations of the system depicted are contemplated. 
For example, the algorithm can be applied to pulsed laser 
ring-down systems, as well. The system described is 
provided only as an example and not as a limitation with 
respect to the ring-down system contemplated by the present 

15 invention. 

Turning to Figure 1, an example of a prior art cavity 
ring-down spectroscopy (CRDS) system is shown. Light is 
generated from a narrow band, tunable, continuous wave diode 
20 laser 20. Single mode diode lasers are well known in the art 
and are produced and commercialized for use in instruments 
such as described herein. The laser 20 can be a single mode, 
continuous wave diode laser tunable in a particular 
wavelength region. 

25 

Laser 20 may alternately be one of a number of 
different types of lasers such as distributed feedback (DFB) 
laser, an "external cavity" diode laser or an optical fiber 
laser. Laser 20 can be temperature tuned by a temperature 
30 controller 30. 

An isolator 40 is positioned in front of and in line 
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with the radiation emitted from laser 20. Isolator 40 
provides a one-way transmission path, allowing radiation to 
travel away from laser 20 but preventing radiation from 
traveling in the opposite direction. Thus, isolator 40 
5 protects laser 20 from back reflections or optical feedback, 
which tend to increase laser noise. 

The light emitted from laser 20 must be coupled as 
efficiently as possible into the optical fiber 48. A single 

10 mode fiber coupler ("F.C.") 50 is provided for that purpose. 
Fiber coupler 50 is positioned in front of and in line with 
isolator 40. Fiber coupler 50 receives and holds optical 
fiber 48 and directs the radiation emitted from laser 20 
toward and through a first lens 46. First lens 46 collects 

15 and focuses the radiation. Because the beam pattern emitted 
by laser 20 does not perfectly match the pattern of light 
propagating in optical fiber 48, there is an inevitable 
mismatch loss. Fiber coupler 50 reduces this loss to about 
3db. 

20 

The laser radiation is approximately mode-matched into 
a ring-down cavity ("RDC") cell 60. A reflective mirror 52 
directs the radiation toward a beam splitter 54. Beam 
splitter 54 directs at least 50%, and typically about 90%, 

25 of the radiation through a second lens 56. Second lens 56 
collects and focuses the radiation into cell 60. The 
remaining radiation, typically 10%, passes through beam 
splitter 54 and is directed by a reflective mirror 58 into a 
water calibration cell 90- -also referred to as a wavelength 

30 select meter. 

Calibration cell 90 is used to tune laser 20. The 



Page 11 of 28 



radiation, which is transmitted through calibration cell 90 
is directed toward and through a fourth lens 92. Fourth lens 
92 is aligned between calibration cell 90 and a second 
photodetector 94 (PD 2) . Photodetector 94 provides input to 
5 computer and control electronics 100. 

Cell 60 is made from two, highly reflective mirrors 62, 
64, with radii of curvature of about 1 meter and separated 
by a distance, d, of about 1 meter. This distance was chosen 
10 based upon stock available mirrors; shorter cells could 
easily be built for a f ield-deployable system with only a 
modest decrease in sensitivity. 

The length of cell 60, L, is a compromise between 
15 opposing requirements. Because RDC decay time is dominated 
by mirror losses and transmissivity, the cavity length 
should be maximized. On the other hand, the maximum length 
is determined by space considerations and by the requirement 
that the cavity be stable with respect to beam size. 
20 According to the theory of stable resonators, in the case of 
two mirrors with the same radius of curvature, r, the 
transverse modes of the resonator have a finite size if the 
cavity length is less than 2r (their size diverges for this 
value) . 

25 

Mirrors 62, 64 constitute the input and output windows 
of cell 60. The sample gas under study flows through a 
narrow tube 66 (1 cm 2 cross section) that is coaxial with 
the optical axis, a, of cell 60. Mirrors 62, 64 are placed 
30 on adjustable flanges or mounts (not shown) that are sealed 
with vacuum tight bellows (also not shown) to allow 
adjustment of the optical alignment of cell 60. 



Page 12 of 28 



Mirrors 62, 64 have a high-reflectivity dielectric 
coating and are oriented with the coating facing inside the 
cavity formed by cell 60. A small fraction of laser light 
5 enters cell 60 through front mirror 62 and "rings" back and 
forth inside the cavity of cell 60. The mechanical alignment 
of the mirrors is most forgiving when the RDC is close to 
confocal (L=r) and becomes extremely critical when either 
the concentric (L=2r) or planar (L is much less than r) 
10 limits are approached. 

Light transmitted through rear mirror 64 (the 
reflector) of cell 60 is directed toward and through a third 
lens 68 and, in turn, imaged onto a first photodetector 70 
15 (PD 1) . Photodetectors 70, 94 may be thermoelectrically 
cooled. 

As shown in Figure 2, an exemplary embodiment of the 
present invention can be outlined as follows: 201: Ring-down 

20 data are digitized by an A-D converter and read into a 
computer as a one -dimensional array. 202: The last hundred 
elements or so are averaged to obtain a DC offset and an 
estimate of the RMS noise. 203: The baseline offset is 
subtracted from the initial data to form a new array. 204: 

25 Starting at the first elements of the new array the data is 
tested for a zero value or a sign change. 205: The last 
point before such a change is the maximum array element that 
can be used in the logarithmic transformation. 206: The 
maximum array element number is scaled by a factor greater 

30 than unity, and array elements greater than this value are 
fit to a straight line and a new estimate of the RMS noise 
determined. 207: The line is extrapolated to the beginning 
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10 



of the record and subtracted from the initial data to 
produce a new array. 208: Once again the array is examined 
for a sign change to obtain the maximum abscissa for the 
logarithmic transformation. 209: After removal of the 
initial few points in the switching transient, a weighted 
least squares fit of the logarithmic transformed data to a 
straight line determines the decay time and the amplitude of 
the ring-down signal and provides estimates of the 
precision. 

Steps 2 06-208 provide the improved performance in the 
presence of low frequency noise as taught by the present 
invention. 

15 As an example of the present invention, a model of a 

ring-down signal with noise is created to mimic experimental 
data obtained from a photodetector and preamplifier applied 
to an A to D converter. It is then fit with the novel 
algorithm of the present invention that subtracts a line 

20 that is fit to the background before applying the standard 
algorithm and compares that fit with one that simply uses a 
constant for baseline subtraction as taught in the prior 
art. Exponential decays of tau = 80 microseconds are used, a 
peak signal is chosen, Gaussian noise is added, a dc offset 

25 is added and then sixty Hz noise is added. The data is 
sampled every microsecond. Lastly with this data, one 
iteration of the prior art fit is used and the present 
invention is used and then the results are compared. 

30 To produce the model data, first, an arbitrary integer 

is selected for the pseudo random number generator seed. 
Next, values for the Gaussian white noise in millivolts, the 
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signal amplitude, the 60 -Hertz RMS level, and the DC offset 
voltage are provided along with the decay time (tau) . In 
one embodiment the inverse of the decay rate, a decay time 
of 80 microseconds is used. The 60-Hertz RMS level must 
5 also be converted to amplitude from an RMS value. 

Next the records are fit to the exponential, first 
using a weighted log transformation. To do this, the 
average value for the last 100 points or so of each decay is 

10 determined to establish a baseline. For these points, the 
signal has fallen off many orders of magnitude below the 
noise. Once the baseline is established, it is subtracted 
from the data and sigma is determined from those last 100 or 
so points, where sigma is calculated from the baseline 

15 fluctuation. The first non-positive array element is 
identified to determine the range of data that can be 
transformed logarithmically and to first order the time for 
the exponential decay to reach the noise level . 

20 Now, data for determining the line that is to be 

subtracted from the initial ring-down data is selected from 
those data points beginning at twice the time that the 
exponential decays to the noise level and continuing to the 
end of the record. This data begins at an array number 

25 designates as (i 0 ) . In a preferred embodiment, no less than 
2 00 points should be used for the fit, which may reduce the 
factor of two in time by which the data is selected. At this 
point, the ring-down signal is as far below the noise as it 
was above the noise initially. This data is fit to a 

30 straight line, which is then subtracted from the initial 
data rather than just a constant as is taught in the prior 
art. The mathematical formula for a straight line is y = 
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(ra)x + b, where y is the ordinate and x the abscissa, and 
the parameters m and b are the slope and intercept 
respectively. In the present case, the ordinate is the array 
of points, s(i), and the abscissa the array number, (i) , 
5 beginning at the value (i 0 ) and continuing to the end of the 
record. 

Once the line is determined, it is subtracted from the 
initial data. This includes finding the new maximum value 

10 of i (determined by the first non positive point) the 
running variable of the digitized record. Next, the 
logarithm of the reduced data is found up to the maximum 
value of i. This transformed data is used to calculate the 
slope and intercept of the line with a weighted average of 

15 the points, as described by Bevington. In addition, the 
reduced average chi 2 and sigma tau are calculated. Lastly, 
the decay rate is to be converted to decay time or tau. 

At this point, new data to produce a new ring-down 
20 curve for an ensemble is provided using the same procedure 
as before with of course new and independent values for the 
noise. Appendix B displays sequentially the calculated 
results for each element of the ensemble of ring-downs; the 
time constant, tau; the amplitude, amp; the number of points 
25 used in the fit, num points; the standard deviation of the 
noise, sigma; the reduced chisq; the standard deviation of 
the time constant, sigma tau; and the standard deviation of 
the amplitude, sigma amp. 

30 Once these values are all determined, this process can 

be repeated using a constant baseline for subtraction 
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instead of the method of the present invention which 
establishes a line to fit for baseline subtraction. 

Given the parameters for the dc offset, decay time, 
5 signal amplitude, 60 -hz amplitude, RMS random Gaussian 
noise, and pseudo random number seed, an array, s(i), that 
represents the ring-down voltage is returned. As can be 
seen in the results of Appendix B, comparing the results 
with the weighted best fit line values for baseline 
10 subtraction to the values generated using a constant 
baseline subtraction to determine decay rates, the present 
invention provides a marked improvement over the prior art. 

Figures 3-12 depict a ring-down data signal sequence 
15 data model. For each figure the y-axis of the depicted data 
models represents signal voltage in arbitrary units and the 
x-axis represents time in units of seconds. 

Figure 3 depicts a measured ring-down signal for the 
20 present data model . The measured signal includes the sum of 
the ring-down signal sequence, a 60 -Hz sinusoidal noise 
signal and broadband random noise signal. Here, it can be 
seen that the total time for the sum of the ring-down signal 
sequences is about 0.016 sec, which corresponds to the full 
25 period of the sinusoidal noise signal, while each individual 
ring-down sequence period is substantially shorter. Each of 
the component parts is depicted separately in Figure 4, 
while the sum of each individual ring-down sequence is 
depicted in Figures 5-7. 

30 

Turning now to Figure 4 there is shown an exemplary 
data model in its component parts including a ring-down 
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signal sequence 401 contaminated with both 60 Hz noise 402 
and broadband random noise 403. As measured in arbitrary 
units, ring-down signal sequence 401 has an -amplitude of 
approximately 100 units with a decay time of 80 microseconds 
5 and a 30 unit DC offset. The SO Hz noise signal has an 
amplitude of 5 units and a DC offset of 15 units. The 
broadband random noise signal has a unit standard deviation 
of 1 . As can be seen in the example depicted in Figure 4, 
the period of the 60 Hz noise is substantially longer than 
10 that of the decay time of a ring-down signal . 

Each of the separate ring-down sequence signals 
depicted in Figure 3 is shown in expanded view in Figures 5, 
6 and 7 . Figure 5 depicts a typical record for 1 ms on the 

15 ascending slope of the sinusoidal noise signal of the first 
ring-down sequence signal for the time interval of about 
0.001s to 0.002s. Figure 6 depicts a typical record for 1 
ms near the peak of the sinusoidal noise signal of the 
second ring-down sequence signal for the time interval of 

20 about 0.004s to 0.005s. Figure 7 depicts a typical record 
for 1 ms on the descending slope of the sinusoidal noise 
signal of the third ring-down sequence signal for the time 
interval of about 0.01s to 0.011s. 

25 It is shown in Figures 5-7, that the algorithm of the 

present invention can correct for low frequency noise 
regardless of the exact frequency of that noise or even if 
it is made up of several frequency components. That is to 
say that Figures 5-7 show that once the ring -down signal has 

30 reached a level where the signal to noise ratio approaches 
values less than 1:1, meaning that only the noise remains, 
the data points can be measured at any point along the 



Page 18 of 28 



sinusoidal noise signal whether it be ascending, descending 
or at its peak value. This can be done because the overall 
slope of the record at any point along the sinusoid 
approximates a straight line for the short time frame in 

5 relation to the entire time for a full period. Accordingly, 
whether data is extracted on the ascending, descending, or 
peak of the sinusoidal noise signal the ability to obtain 
more precise time decay results is not affected. Because 
the time segment of the record is so small the addition of 

10 the sine signal section closely approximates a straight 
line. Therefore, one can use a straight line approximation 
beginning from the point of i 0 to the end of the record to 
subtract out the low frequency noise. 

15 Turning now to Figure 8, there is again shown the first 

ring-down signal in expanded view. In this figure, the last 
one hundred points of the record are used to find the 
average DC offset to subtract from the ring-down data. To 
find the average DC offset, the average value of the final 

20 100 data points of the measured signal is calculated. The 
final 100 data points are used, because the ring-down signal 
will have reached a minimum value at this point. Once the 
average DC offset is calculated that value is subtracted 
from the measured signal data points, thereby shifting the 

25 measured signal along the y-axis. 

Turning now to Figure 9, there is shown the ring -down 
signal of Figure 8 after the DC offset has been subtracted. 
The points from the beginning of the shifted record to the 
30 point 901 before the first zero or negative value is used in 
the known algorithm to determine the exponential decay rate 
or time. According to that method, in this example the last 
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data point that can be used in the logarithmic 
transformation 901, which corresponds to the point before 
the first zero or negative point is the 276 th point of the 
record. Here, the last hundred points of the record has an 
5 average value of 4 8.491, which is the DC offset to subtract 
from the data. 

Figure 9 provides an expanded view of Figure 5, after 
the DC baseline has been subtracted. The points from the 

10 beginning of the record to the point before the first zero 
or negative value, are used in the standard algorithm to 
determine the exponential decay rate or time. Here, the 
first negative value is recorded at 0.001277 seconds, or at 
the 277 th point of the record. The time value for the 

15 positive point of the record establishing the point before 
the first zero or negative value is doubled 902. In this 
example, that time value, 0.001276 sec. is doubled to 
0.001552 point of the record, which corresponds to the 552 nd 
point, 1001, as shown in Figure 10. Once that point is 

20 determined a straight line is fit to the end of the record 
for that value 1001. In Figure 11, that value, 1001 which 
corresponds to 1101, is then used to find a line 1102 that 
will be subtracted at each point from the original data as 
shown in Figure 11. In Figure 12, the ringdown signal of 

25 Figure 11 is shown in an expanded view after the fitted line 
has been subtracted from it. The last positive value is 
again determined, this time for use in a logarithmic 
transformation. In this example, the last positive point is 
at 326. Figure 13 provides a data model of Figure 12 for 

30 the first 326 points, up to the last positive value, plotted 
on a logarithmic scale. Figure 13 also shows how a straight 
line fit would determine the exponential decay rate of the 
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ringdown signal . The points are to be weighted in 
proportion to their amplitude in the fitting process, 
accentuating the higher valued points. The techniques of 
the present invention described herein may be implemented in 
5 hardware or software, or a combination of the two, but are 
not limited to any particular hardware or software 
configuration. The foregoing descriptions of the specific 
embodiments of the present invention have been presented for 
purposes of illustration and description only. They are not 

10 intended to be exhaustive or to limit the invention to the 
precise forms disclosed, and obviously many modifications 
and variations are possible in light of the above teaching. 
The embodiments were chosen and described to best explain 
the principles of the invention and its practical 

15 application, to thereby enable one of ordinary skill in the 
art to best utilize various embodiments with various 
modifications as are suited to the particular use 
contemplated. It is intended that the scope of the 
invention be defined by the claims appended hereto and their 

20 equivalents. 
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